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The tree-level two-point amplitudes for the transitions jf j'f'•, where / is a fermion 
and j is a generalized current, in a constant uniform magnetic field of an arbitrary 
strength and in charged fermion plasma, for the jff interaction vertices of the scalar, 
pseudoscalar, vector and axial-vector types have been calculated. The generalized current 
j could mean the field operator of a boson, or a current consisting of fermions, e.g. the 
neutrino current. The particular cases of a very strong magnetic field, and of the coherent 
scattering off the real fermions without change of their states (the “forward” scattering) 
have been analysed. The contribution of the neutrino photoproduction process, ye 
euP, to the neutrino emissivity has been calculated with taking account of a possible 
resonance on the virtual electron. 

Keywords: Charged fermion plasma; magnetic field; Landau levels; astrophysics. 
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1. Introduction 

Nowadays, there exists rather keen interest to astrophysical objects with the scale of 
the magnetic field strength near the critical value of B^. = m^je ~ 4.41 x 10^^ G[^ 
This group of objects includes the radio pulsars and the so-called magnetars, which 
are the neutron stars featuring the magnetic field strengths from 10^^ G (radio pul¬ 
sars) to 4x 10^^ G (magnetars), see Ref.[T]and the papers cited therein. The spectra 
analysis of these objects also provides an evidence for the presence of electron- 
positron plasma in the radio pulsar and magnetar environment, with the minimum 
magnetospheric plasma density being of the order of the Goldreich-Julian densityP 



( 1 ) 


®'We use natural units c = h = k^ = 1, rUe is the electron mass, and e is the elementary charge, 
mf and Cf are the fermion mass and the fermion charge. 


1 







January 15, 2015 1:28 WSPC/INSTRUCTION FILE 

Kuzn’ Rum' Shlen' Gen' Compt' corr 


2 A. V. Kuznetsov, D.A. Rumyantsev, D.M. Shlenev 

where P is the rotational period. It is well-known that strong magnetic field and/or 
plasma could have an essential influence on various quantum processes,because 
the external active medium catalyses the processes, by changing their kinematics 
and inducing new interactions. Therefore, the effects of magnetized plasma on mi¬ 
croscopic physics should be incorporated in the magnetosphere models of strongly 
magnetized neutron stars. In the present paper we consider the two-point processes, 
because such reactions can have possible resonant behavior, and therefore they could 
be very interesting for astrophysical applications.^ 

The investigation of the two-point processes in an external active medium (elec¬ 
tromagnetic field and/or plasma) has a rather long history. The most general expres¬ 
sion for a two-vertex loop amplitude of the form j —>■//—>■/' in a pure constant 
uniform magnetic field and in a crossed field was obtained previously in Ref. |8l 
where all possible combinations of scalar, pseudoscalar, vector, and axial-vector 
interactions of the generalized currents j and j' with fermions were considered. 
The generalized current j could mean the field operator of a boson, or a current 
consisting of fermions, e.g. the neutrino current. 

The typical example of a tree-level process with two vector vertices in the pres¬ 
ence of magnetized plasma is the Compton scattering, ye —>■ ye, as a possible channel 
of the radiation spectra formation. In this case, both generalized currents j and j ' 
mean the photon field operators. This process was studied in a number of papers, 
see e.g. Refs. IMH but the results were presented there in the form without taking 
account of the photon dispersion properties. In Ref. [121 this neglect was corrected. 
The expression for the Compton scattering amplitude, with the initial and final 
electrons being on the lowest Landau level was presented in Ref. [12] in the explicit 
Lorentz invariant form. The other example of the Compton like process with the 
vector and axial-vector vertices, the photon transition into the neutrino pair in the 
presence of magnetized plasma, ye evv, was studied in Ref. |T71 In this case, in 
our terms, the initial generalized current j means the photon field operator, while 
the final generalized current j' means the neutrino current. The local limit of the 
weak interaction is supposed to be valid. One of our goals in this paper is to improve 
the approach of Ref. [TTl in order to present the results in a manifestly covariant 
form. Additionally, as we believe, our results would be better applicable for an anal¬ 
ysis of the other photon-fermion scattering processes with the production of exotic 
particles, such as axion, neutralino, etc. 

Thus, we consider the tree-level two-point amplitude for the transition of the 
type jf —>■ j'f with the intermediate virtual fermion state. The analysis is per¬ 
formed in a constant uniform magnetic field and charged fermion plasma, for differ¬ 
ent combinations of the vertices that were used in Ref. |8l Particularly, we generalize 
the results obtained in Ref. |H1 to the case of magnetized plasma, since such a sit¬ 
uation looks the most realistic for astrophysical objects. Such a generalization was 
performed in part in Ref. [T8| for the case of the photon polarization operator in a 
magnetized electron-positron plasma. 






January 15, 2015 1:28 WSPC/INSTRUCTION FILE 

Kuzn’ Rum' Shlen' Gen' Compt' corr 


Generalized two-point tree-level amplitude ... 3 


The paper is organized as follows. In Sec. we calculate the scattering ampli¬ 
tudes for different spin states of the initial and final fermions. We present here only 
the amplitudes for the case when both vertices are of the pseudoscalar type. The to¬ 
tal set of the amplitudes for the jff interaction vertices of the scalar, pseudoscalar, 
vector and axial-vector types, in a constant uniform magnetic field of an arbitrary 
strength and in charged fermion plasma can be found in the extended paper.l^ All 
the amplitudes are presented in the explicit Lorentz and gauge invariant forms. The 
application of the obtained results to the calculation of the neutrino photoproduc¬ 
tion process amplitude and other characteristics in the resonant case is given in 
Sec. 1^ Final comments and discussion of the obtained results and possible astro- 
physical applications are given in Sec. In [Appendix A[ we present the fermion 
wave functions used in our analysis, namely, the solutions of the Dirac equation in 
external magnetic field, being the eigenfunctions of the magnetic moment operator. 
In the next two Appendices, we present the expressions for the amplitudes in the 
special cases where they can be essentially simplified. In [Appendix B[ we consider 
the particular case, when the initial and final fermions occupy the ground Landau 
level (the strong field limit), for all types of the jff interaction vertices. A coherent 
scattering of neutral particles off the real fermions without change of their states 


(the “forward” scattering) is analysed in Appendix C 


2. The set of expressions for the amplitudes 

The generalized amplitude of the transition jf —>■ j'f' will be analyzed by using 
the effective Lagrangian for the interaction of a generalized current j with fermions 
in the form 


C{X) = Y,gk[^fiX)Tk^f{X)]MX), (2) 

k 

where the generalized index k = A numbers the matrices P/j.: r 5 = 1, Pp = 

75, Fy = 7a, Ta = 7 a75; jfe(A) are the generalized currents (j's, jp, jva or Jau) or 
the field operators of single particles, e.g. of the photon or axion, see below, gp are 
the corresponding coupling constants, and 4'y(X) are the fermion wave functions. 

Indeed, using the Lagrangian ([^, one can describe a large class of interactions. 
For example, it may be: 

i) the Lagrangian of the electromagnetic interaction, when k = V, gv = ~e/, 
Fyjy = is the four-potential of the quantized electromagnetic field: 

C{X) = -ef [^f{X)j^A^{X)^fiX)] ; (3) 


ii) the Lagrangian of the fermion-axion interaction, when k = A, gA = Cf/ (2 fa), 
^aJa = l^l5d^a{X), a{X) is the quantized axion field, fa is the Peccei-Quinn 
symmetry violation scale, Cf is the model dependent factor of order unity: 


^X) = P[^fiX)rj,^fiX)]d^a{X), 

^Ja 


( 4 ) 
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Fig. 1. The Feynman diagrams for the reaction jf j' f' • Double lines mean that the effects of 
an external field on the initial and final fermion states and on the fermion propagator are exactly 
taken into account. 


iii) the effective local Lagrangian of the four-fermion weak interaction, when 
k = Y^ gy = G^Cy!^/2 and k = qa = ~GyGa! 

C{X) = ^ [§/(X)7„(Cy - CAlB)^f{X)] UX), (5) 

where Ja{X) = j/(X) 7 „(l — 'y^)i'{X) is the current of left-handed neutrinos; Gy = 
±1/2 ± 2sin^ 9w, Ca = ±1/2, and 9y/ is the Weinberg angle. Here, the upper sign 
corresponds to neutrinos of the same flavor f {v = Vf), when there is an exchange 
reaction both of W and Z bosons. The lower sign corresponds to the case of another 
neutrino flavors (z/ ^ Vf), when there is only Z boson exchange. The conditions of 
applicability of the effective Lagrangian © should be specified. First, it is the 
condition of relatively small momentum transfers, \q'^\ <C , where mw is the 
W boson mass. And second, the condition that additionally arises in an external 
magnetic field, is eB <C . We will consider physical situations where both of 
these conditions are satisfied. 

In a general case with the Lagrangian ([^, the S'-matrix element in the tree 
approximation is described by the Feynman diagrams shown in Fig. [^and has the 
form 


QS S _ 

' k 


gkge J d^Xd^Y {jk{X)jk>{Y)) 

(jfc, Tfc <->■ j'j., F/). 


^!;„(y)Ffc,5(y,x)r,vE-w(A) (6) 


Here, = {Ei,p) and p'^ = (if/, ,p') are the four-momenta of the initial and 
final fermions correspondingly, X^ = {Xq, X 1 A 2 , Xs), the currents between the 
angle brackets mean the matrix element between the corresponding initial and final 
states, d/p i{X) are the fermion wave functions in the presence of external magnetic 
field, where the subscript describes a state with definite components of the four- 
momentum p and with the Landau level number i, while the superscript describes 
the spin state s. 

There exist several descriptions of the procedure of obtaining the fermion wave 
functions in the presence of an external magnetic field by solving the Dirac equation, 
see e.g. Refs. I^DH^ and also Refs. 00 In the most cases, the solutions are presented 
in the form with the upper two components of the bispinor corresponding to the 
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fermion states with the spin projections 1/2 and -1/2 on the magnetic field direction. 
Here, we have found it more convenient to use another representation of the fermion 
wave functions, being the eigenstates of the magnetic moment operatorl2£lE^ Some 
details on these wave functions are presented in Appendix A[ 

The currents jp in Eq. ([^ can be expressed through the amplitudes in the 
momentum space: 


JkiX) = 


g-i(gX) 

7 ^ 


ik{q) ■ 


(7) 


We use the fermion propagator in the form of the sum over the Landau levelslUlZl 


OO 


S{X,X') = ^ SniX,X'), 

n—0 


( 8 ) 


^ r dpodpydp^: e ^^^11 

J (27r)3 — TTij — 2/3n-I-ie 

X exp |- ^ -Pj, [ Ai + A/ - i (A 2 - A 2 ')]| 
X I [(p 7 )„ +mf] [i 7 _ iJn(e') -b n+ 2 n 


-bi 


2n 7^77- 1 (0 (e') - 7T+ (e) iJn-1 (C')] I, 


(9) 


where ^ is defined by Eq. (A.14) and is obtained from ^ by substituting Xi —>• X(. 


Hereafter we use the following notations: four-vectors with the indices J_ and 
II belong to the Euclidean {1, 2} subspace and the Minkowski {0, 3} subspace 
correspondingly. Then for arbitrary 4-vectors A^, one has 


= (0, Ai, A 2 , 0), A/ = (Ao, 0,0, A 3 ), 

(AB)^ = (AAB) = AiBi -k A 2 B 2 , 

(AB)ii = (AAB) = A 0 B 0 -A 3 B 3 , 

where the matrices A^^, = A^,^ = are constructed with the di¬ 
mensionless tensor of the external magnetic field, (p^,^ = and the dual 

tensor, ■ The matrices A^,^ and A^,^ are connected by the relation 

AA“^ — = diag(l, —1, —1, —1), and play the roles of the metric tensors in 

the perpendicular (T) and the parallel (||) subspaces respectively. 
















January 15, 2015 1:28 WSPC/INSTRUCTION FILE 

Kuzn’ Rum' Shlen' Gen' Compt' corr 


6 A. V. Kuznetsov, D.A. Rumyantsev, D.M. Shlenev 


After integration in Eq. (|^ over di^X and d'^Y we obtain 

i(27r)3j(3)(P-p'-g') 


QS S _ 

^k'k — 


Mi,I 


(10) 


^2qoV2q^V2EeLyL,2E',LyL, 

where 5^{P - p' - q') = 5{Po - E', - q^)S{Py - Py - qy)5{P^ - p'z - 9z')> Pa = 
{p + q)a, a = 0,2,3, and the partial amplitudes Mf.,% can be presented in the 
following form: 

— exp [—id] 


‘fc'fc 


2 ^MiMi,{Mi + mf){M£, +'rnf) 


( 11 ) 


jexp 

'iiqpq'Y 

Qy + Ida; 

-i 

1- 

1 

1 _ 

[ 2/3 J 

. nA? . 


-1 

_ 1 


E 

n—O 


{qAq') - i{q(pq') 


yjp (1)S 'S 
'^k'k 




2^/2 


P|^ — nij — 2f3n 


+(-1)^+^ exp 


X V f (g^gQ + l(g‘^g') 
Po V ^ 


iiqpq'Y 

1 

+ 

1 _ 

-t 

1 - 

1 

1 _ 

[ 2/3 J 

\A? _ 


. . 


-i’ 


n 


(2)s 's 


kk' 




P|i'2 — ruj — 2f3n J ’ 


where 0 = {q^c — q'x){py +Py)/(2/3) is the general phase for both diagrams in Fig. 

P'a = {p-q')a- 

The main part of the problem is to calculate the values * which are 

expressed via the following Lorentz covariants in the {0, 3}-subspace 


{pAp') + MiMi 


{pAp') + MiMt 


- |M^(Ap')a + M^/(Ap)a| , 
{Mt{ipp')a+ Mt,{pp)a} , 


ICia — 


fdla — 


/C3 = -y/2 {pAp')+MiMi, 

/C 4 = -, / —- {ppp') ■ 

V (pApO + ^ 

The following integrals appear in the calculations: 

1 


( 12 ) 

(13) 

(14) 

(15) 


^ / dZe-^ Hr 

/TT ,/ 






Qy + igx 
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where, for n ^ £ 

TnA^) = = (-l)""^ J„y(cr), (17) 

and L^{x) are the generalized Laguerre polynomials.^ 

Below, the results are presented for the values * in the case when both 

vertices are of the pseudoscalar type, k = k' = P. The total set of the values 
jff interaction vertices of the scalar, pseudoscalar, vector and 
axial-vector types, in a constant uniform magnetic field of an arbitrary strength 
and in charged fermion plasma is presented in the extended paper.^ 

Hereafter we use the following definitions: = In,e {q±/{‘^f3)) and = 

^n,t' {q±/{‘^P))- For definiteness, we further consider the fermion with a negative 
charge, e/ = — |e/|. 

In the case when j and j' are the pseudoscalar currents (k = k' = P) we obtain 
= -5p5^JPjp|2/3v^[(-*CiP) 

+ {Mi + ruf) [[ICiP) - TO//C 3 ] (18) 

-2/3Vn/C3 [v^(M^/ -I- P \/V{Mi + | ; 

'R!'pp~ = -igp5pJpjp|\/2/^£' {Mi +mf) [{IC 2 P) + mflC4]l^i,In,e 
-,/We {Mp + ruf) [{IC 2 P) - mflC4Wu'Wn-i,t-i (19) 

-,/Wn lC/{Mi +mf){Mp + m/lW/'-W/ - 2/3v^X^y/X„_i,f_i] |; 

l/pp^ = igpg'pjpjp^y/We {Mp P ruf) [(/C 2 P) - mfWl/Wmf- 
-Wm' {Mi P iTif) [{IC 2 P) P mfWlL-ix'Wn-i,i-i (20) 

-,/Wn lCi[2p,/lrXW,pWn,i - {Ml P mf){Mi> P mW/pX^-i^i-i] | ; 
Pp]r“ = -gpgpjpjp^{Mi Pmf){Mp P nif) [{JCiP) - mfWX/Wn,^ 
P2fi/Jl' [{ICiP) P mfIC3]X^_j^^^,_^Xn-i^i-i 

—2l3/nlC3 ['s/G {Ml P mf)X^_^i,_^Xn^i P '/e{Mp P mf)X/i,Xn_ii_ij . 


(21) 








January 15, 2015 1:28 WSPC/INSTRUCTION FILE 

Kuzn’ Rum' Shlen' Gen' Compt' corr 


8 A. V. Kuznetsov, D.A. Rumyantsev, D.M. Shlenev 


To obtain the contributions T^pp * from the second diagram of Fig. |lj the 
following replacements should be made in Eqs. (181—(21): Pc 


P' T 4-^7' 

_ __ ^^ a; -^rn,n 

The results obtained for the case of the magnetic field of an arbitrary strength 
can be essentially simplified in several special cases. In [Appendix B[ the set of 
expressions for the amplitudes in the limit of relatively strong field is presented, 
where the initial and final fermions are on the ground Landau level, £, £' = 0, but 
the virtual electron can occupy an arbitrary Landau level, n ^ 0. 

One more case when the amplitudes can be essentially simplified is the process 
of a coherent scattering of the generalized current j off the real fermions of magne¬ 
tized plasma without change of their states (the “forward” scattering). The set of 
expressions for the amplitudes in this case is presented in [Appendix C| 


3. Neutrino luminosity 

As an illustration of the results obtained, let us construct the amplitude of the 
neutrino-antineutrino pair photoproduction, ye -7 euv, in a strongly magnetized 
cold plasma when the temperature T is the smallest parameter of the problem, i.e., 
T <C /Te — vrie {pe is the chemical potential of the electron gas, rrie is the electron 
mass) with taking account of a possible resonance on a virtual electron. 

At the same time, our main goal is to obtain the expression for the neutrino 
emissivity caused by the process ye -7 evv. In turn, the neutrino emissivity can be 
defined as the zero component of the four-vector of the energy-momentum carried 
away by the neutrino pair due to this process from a unit volume of plasma per 
unit time. Here, we neglect the inverse effect of the energy and momentum loss on 
the state of plasma. The neutrino emissivity can be represented in the formP^ 


Q 




-T. ^ 


d?k 


(27r)3 2qo 


^(fo) 


d^p 


( 2^)2 2Ei 


fe{Ee) 


X 


h _ f d^p2 , 

{2nY2E^, ^ (27r)3 2£i {2n)^2e2^° 


X (27r)^ 6^{P -p' - q')\M-,,e^enA'^ , 


( 22 ) 


where f-y{qo) = — l] is the equilibrium distribution function of an ini¬ 

tial photon with the four-vector q>^ = (go,k); fe{Ee) and fe{E^,) are the equi¬ 
librium distribution functions of initial and final electrons, respectively, fe{Ei) = 

_|_ £2 is the neutrino pair energy, £ 1^2 = |Pi, 2 |; 
d^p = dpydpz; V = LccLyLz. is the plasma volume. 

In calculating the amplitude of the process ye -7 evv, we consider 

the case of relatively small momentum transfers compared with the W boson mass, 
|g'2| ^ . Then the corresponding interaction Lagrangian can be written as 
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follows, see Eq. ([^: 


£ = ^ [5'7q(Cv - C'a 75)^'] [i^7a(l - 75 ) 1 ^] + e(^'7Q5') , (23) 

where A^, is the four-potential of the photon field. 

Comparing (23) with the Lagrangian of the general form (§ we find that the 
amplitude squared of the process qe —)■ eiyp can be represented as: 


\M 


'ye—^eui' \ 




(24) 


and in formulas ( |B.2 1, ( |B.17[ )" p3.20[ ) one should put mj = rrie, gy = e, where 
e > 0 is the elementary charge, ja = £a is the initial photon polarization vector, 
9v = GfCvIV^, g'A = -G-pCaIV^, ja = i^7a(l - 75)i^- 

It should be noted that the virtual electron resonance occurs only in the s channel 
diagram (the first diagram in Fig. [^. Nevertheless, even with the simplification 
caused by the resonance behavior, the problem under consideration is still enough 
cumbersome, because the charged fermions can occupy arbitrary Landau levels. The 
problem could be significantly simplified in the physical conditions of magnetars. 
Indeed, in the outer crust of a magnetar, the following hierarchy of parameters 
should existP^ eB T^. Thus, the electron plasma can be considered as 

a strongly magnetized one, and under these assumptions one can approximately 
assume that the initial and the final electrons would occupy the ground Landau 
level (£ = £ = 0), while the virtual electron can occupy an arbitrary Landau level. 


In our case, s' = s = — 1 and the amplitude squared (24) takes the form 


\M 


-ye—^euu \ 


E 

n—l 


{qAq') - ijqipq') 


TZr, 


— ml — 2j3n 


(25) 


where 


1 

nl 


^2 N 

"/2 /n 


I ^ 

2/3, 

) 

>(i) 

'^vv 

and Pav 


1/2 


exp 


qi + q'r 


4/3 


p(l) I p' 
liyy + -it 


( 1 ) 

AV 


(26) 


To accurately take into account the resonance behavior in the process ye —t evv, 
it is necessary to calculate radiative corrections to the electron mass, caused by the 
combined action of a magnetic field and plasma. This calculation is a separate chal¬ 
lenge. However, because of the smallness of these corrections, we can approximately 
replace ml —>■ ml — iPorn in the denominator of Eq. (25). 


As it was already noted, the main contribution to the amplitude arises from the 
resonance region, so that we can approximately replace the corresponding part of 
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Eq. (25) by the d function: 


OO 

n—1 


|7^nP 

(p2-m2-2/3n)2 + p2r2 


-E 


PnF 


n 




2pn)\T^nW 


(27) 


where r„ is the total width of the change of the electron state. This width can be 
represented in the forrrP^ 


r — r 

■L n — -L 


abs 


■ r 


eo7— 




(28) 


Here 


per 

^ eo7 —>-671 


1 


290(2^)3 2Eo(271)2 


fe{Eo) 


x{27T)H^{P-p")\Meoy^eS 


(29) 


is the width of the electron creation in the nth Landau level. 


With taking account of Eq. (281, the amplitude squared of the process ye 
takes the form: 


evv 


1*^76—I — ^ ) 

n 

d^p 


i2.fSHP-pl 


= E 

n—1 


n—1 

2^// 


(2^)2 2E" 


2i?"r„ 


(271)2 2E" 


UE':) {2^Y5\P-p") 


\Rn 


O fP," ycr 
^-^n ^ eQ-f—>-en 


Here we have used the property of the 6 function: 

J(P|f -m!- 2/3n) = ^ S(Po - P") 


(30) 


(31) 


where P" = + 7n2 + 2/3n. 

On the other hand, in the case of resonance the expression for |P„|2 being 
averaged over the photon polarizations can be factored in the strong field limit 
/3 TOe as follows (see, for example. Ref. ED: 


\Pn? = |JWeo^_,e„ P | AJe„-J.eo!^P P , 


(32) 


where 


I . > |2 Sttg; 

1*'^eo 7 —>-671 I — 6Xp 


n\ 




Ml (pAq) 


xE 


A=1 



2 

+ 



+ a A 


(A) 

3 


2 


(33) 


is the amplitude of the absorption of a photon in the process egy —l e„, when an 
electron passes from the ground Landau level to a higher Landau level n. The pa¬ 
rameter a = {ppq) / {pAq) = ±1 determines the direction of the photon propagation 
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with respect to the magnetic field direction. are the expansion coefficients of 
the photon polarization vector over the basis of the 4-vectors: 

b^p^ = {Aq)p - qp ql, 6^ = qp . (34) 

Given the gauge invariance, one has: 

= (35) 

i-1 

Finally, the amplitude squared of the electron transition from the nth Landau level 
to the ground level with the creation of the neutrino-antineutrino pair takes the 
form: 


\M 




(^2 

|2 _ 

n! 


exp - 




12 


J2 


Ml 


2/3 / \ 2^ / {p'Aq') 


Cvip'Aj') - Ca{p'(P3') 


(j'Aj'*) 


,/2 


Cv{p'Aq')-CA{p'pq')\ - Cv{p'Aq') - Ca{p' m') 
J q'J 


xRe [q'Aj') Cvip'Aj'*) - CAip'pj'l 


(36) 


account Eqs. (29) and we obtain: 


Substituting Eq. (|30|) into the expression for the luminosity (221, and taking into 

(37) 


Q'yeo—^eQi'i/ — ^ ^ Qen-^-eoi^i' ? 
n—1 


where 


1 
Ij'T 


Qen—^eQizo 
d^pi 

^(27r)3 2ei (27r)3 2e2 


/ (27r)2 2EII ^ (27r)2 2E' 


d^P2 


lo (27r)3 J3(p" -p' - q')\M, 


(38) 


is the neutrino luminosity due to the process e„ —> e^vu. This result coincides, up 
to notation, with the result of Ref. [221 


4. Discussion 

In this paper, we have calculated the tree-level two-point amplitudes for the tran¬ 
sitions jf —>■ j'f in a constant uniform magnetic field of an arbitrary strength, 
and in charged fermion plasma, for generalized vertices of the scalar, pseudoscalar, 
vector and axial vector types. It is remarkable, that all the amplitudes obtained are 
manifestly Lorentz invariant, due to the choice of the Dirac equation solutions as 
the eigenfunctions of the covariant operator In this case, partial contributions 
to an amplitude from the channels with different fermion polarization states are 
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calculated separately, by direct multiplication of the bispinors and the Dirac matri¬ 
ces. This approach is an alternative to the method where the amplitudes squared 
are calculated, with summation over the fermion polarization states, and with using 
the fermion density matrices, see, e.g. Refs. ISH ESI However, the use of the density 
matrix in a magnetic field, as is usually done in the absence of a magnetic field, 
in the case of the two-vertex processes leads to extreme difficulties in analytical 
calculations. 

The set of the amplitudes for the transitions J/ —> J in a constant uniform 
magnetic field of an arbitrary strength, and in charged fermion plasma, presented 
in this paper, can be used as a reference book in the investigations of the quantum 
processes in external active media. The field effects are taken into account exactly, 
because exact solutions of the Dirac equation are used. Owing to this, the expres¬ 
sion obtained here for the amplitude is quite general; in particular, it can be widely 
used to analyze various physical phenomena and processes in a magnetic field and in 
plasma. The amplitudes Mss Mpp, which are diagonal in the generalized cur¬ 
rents, differ only in factors from the external-medium-induced contributions to the 
mass operators of the corresponding scalar and pseudoscalar fields. The amplitude 
Mvv defines, for example, the medium-induced part of the photon polarization 
operator. The amplitudes Myv and Mva describe the process amplitude for the 
radiative transition of a massless neutrino p -a v^. Similarly, one can obtain the 
amplitudes for the axion decay a ^ vv and for axion-photon oscillations by means 
of the corresponding substitutions. 

Furthermore, the results obtained can be used to analyze the reactions with a 
possible resonance on the virtual electron (see e.g. Ref. EH)- It is well known that 
the processes of this type play an important role in the magnetospheres of isolated 
neutron stars, providing the production of e+e” plasma.l^ 

Although the obtained formulas look quite cumbrous, there are certain areas of 
their application. We emphasize that these formulas are derived for a general case, 
namely, for arbitrary values of the magnetic field, therefore, they can recover the 
results within the limits of weak and superstrong fields. Further, in this general 
form the formulas definitely may be used for numerical calculations. 
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Appendix A. Solutions of the Dirac equation in an external 
magnetic field 

In this Appendix, we present the fermion wave functions as the solutions of the 
Dirac equation in the presence of an external magnetic field, and simultaneously as 
the eigenfunctions of the magnetic moment operator HUHH 

In Ref. [22l an operator was introduced which was called the generalized spin 
tensor of the third rank. In modern standard notations, the operator takes the 
forn 0 


^ /ii/A 2 (-^A 




(A.l) 


where - 7i^7m)/ 2, and = \d^ - Cf = (ido - Cf Aq, -iV - e/A) 

is the generalized four-momentum operator with A^ being the four-potential of an 


external magnetic field. Taking the component of the operator (A.l) and taking 


into account that in the Schrodinger form of the Dirac equation one has ido = H, 
where = 70 (pfP) A rrif^o + e/Ao is the Dirac Hamiltonian, one can construct 
the vector operator 


/li — 2 ^ijk FjfcO •. 


(A.2) 


where Sijk is the Levi-Civita symbol. This is the magnetic moment operatorj221E^ 
which can be presented in the form 


p = m/S - i7o75[S X P] 


(A.3) 


It is straightforward to show that the components of the operator (A.3) commute 
with the Hamiltonian, i.e. H and pz have common eigenfunctions. In the non- 


relativistic limit, the operator (A.3) is transformed to the ordinary Pauli magnetic 


moment operator, thus having an obvious physical interpretation. 

It appears to be convenient to use the fermion wave functions as the eigenstates 
of the operator 


Pz = TO/Sz - i7o75[S X P]z , (A.4) 

where P = —iV — e/A. We take the frame where the field is directed along the 
z axis, and the Landau gauge where the four-potential is: A^ = (0, 0,xiJ,0). It is 
convenient to use the notation /3 = |e/|iJ, and to introduce the sign of the fermion 
charge as p = e//|e/|. 

Our choice of the Dirac equation solutions as the eigenfunctions of the operator 
Pz is caused by the following arguments. Calculations of the process widths with two 
or more vertices in an external magnetic field by the standard method, including 
the squaring the amplitude with all the Feynman diagrams and with summation 


'^It should be noted that in Ref. 1221 the covariant bilinear forms were constructed of Dirac matrices 
by inserting them not between bispinors 'ip and ip as accepted in modern literature, EU but between 
bispinors ip'^ and ip. 
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or averaging over the fermion polarization states, contain significant computational 
difficulties. In this case, it is convenient to calculate partial contributions to the 
amplitude from the channels with different fermion polarization states and for each 
diagram separately, by direct multiplication of the bispinors and the Dirac matri¬ 
ces. The result, up to a total for both diagrams non-invariant phase, will have an 
explicit Lorentz invariant structure. On the contrary, the amplitudes obtained with 
using the solutions for a fixed direction of the spin, do not have Lorentz invari¬ 
ant structure. Only the amplitude squared, summed over the fermion polarization 
states, is manifestly Lorentz-invariant with respect to a boost along the magnetic 
field direction. 

The fermion wave functions having the form 

„-i(E„Xo-PyX2-p^X3) TTs(f:\ 

vI/« (x) = __ 

+ mf )LyL, ’ 

where 

E„ = y/M^+pl , M„ = \Jrrij + 2j3n , 
are the solutions of the equation 

= s = ±l. 


(A.5) 


(A.6) 


(A.7) 


It is convenient to present the bispinors in the form of decomposition over 

the solutions for negative and positive fermion charge, 


um = 

where 

u--{0 = 


ut-{0 = 


/ -iV2/3npzE„_i(^) ^ 

{Eji + Mn){Mn + mf)Vn{C) 
-iv^(A„ + M„)E„_i(0 

\ -Pz(M„ -f TO/) 14(C) / 

/(A„ + M„)(M„+TO/)y„_i(C)\ 

-V2/3npz; 14(C) 

Pz(M„ -I- to/)14-i(C) 
iv^(A„ + M„)y„(C) / 


(A.8) 


(A.9) 


(A.IO) 


V 
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u-AO = 


ulAO = 


( i^/AnPzyn{0 \ 

{En + Mn){Mn + mf)Vn-l [0 

iAA^{En+Mr,)VA) 

V -p^{Mn+mf)Vn-l{0 ) 

( (F;„ + + m/)K (0 \ 

V2/3np2l4_i(^) 

Pz{Mn + mf)Vn{C) 

V -iv^(E„+M„)K_i(c) y 


(A.ll) 


(A.12) 


(n = 0 , 1 , 2 ,...) are the normalized harmonic oscillator functions, which are 
expressed in terms of the Hermite polynomials 


giiA -en 

va) = ^j=^ha). 

y'2'^n\A 




(A.13) 

(A.14) 


Appendix B. The set of expressions for the amplitudes for ground 
Landau Level, £ — — 0 

In this Appendix, we consider the limit of relatively strong field, where the initial 
and final fermions are on the ground Landau level, £, V = 0, but the virtual 
electron can occupy the arbitrary Landau level, n 7 ^ 0. In this case s = s' = — I, 
Ml = Ml' = ruf, and 


^n,oix) = e ®/ 23 ,n /2 ^ ( 3 ;) = . (B.l) 

Vn! V X 

Denoting = (2 to/)^ R^k''k obtain the following expressions for the am¬ 

plitudes ( 12 ) with the vertices of the scalar, pseudoscalar, vector or axial vector 
types 


Mk'k = -exp[-i 6 »] exp 


/ 2 -| °° 


4/3 


y 1 

77 .! 


X < exp 


' exp 




. 2/3 

i(g'Pg') 

2/3 


n—0 
/\\^ 


( (gAg') -i(gy)g ') 

I 2/3 

(gAgO+KW)^ ” 

2/3 


x^k'k 


P'^ — m'j — 2/3n 
?( 2 ) 


(B.2) 


R 


kk' 


P! 2 — m f — 2/3n \ 
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where 


RU = asg'sjsj's [i^iP) + mfJCz] ; 
R^ss ~ R^ss(a ^ ~a ) > 


R% = gsg'pjsjp [(^2^) - TO//C4] ; 


R^sp = -gsg'pjsjp [(^2-P') + TO//C4] ; 


R^'vl = gsgvjs\ (PAj'O/Ca + (Pw')/C4 + m/(/Cij') 




Rfv = gsgvjs\{P'Rj')K.z - {p'vj'Wa + 

[(^'Aj')+ i(gVj')]A:3|; 


Ras = gsg'AJs\ (pAj')/C4 + {pm^^ - mfilC2j') 


[(^'Aj') - i(gVj')] A :4 


4 a = gsgAJsl iP'^j')IC 3 - (P'Aj')/C 4 - nifilC^j') 


lig'Rj') + Kg'‘pj')]i^ 4 : 

^_L 


(B.IO) 


4 p = -gpg'pjpj'p [i^iP) - PifK-z] ; 


(B.ll) 


4p = 4p(9 ^ - 9 '); 


(B.12) 


4 p = -gpgvjps {PRj')K.A + {P<fj')iC3 + mf{iC2j') 


[(^'Aj') - i((?Vj')] a:4 

n ' ^ 


(B.13) 


4y = gpgvJp\ {P'vj')K.3 - {P'Rj')K.i + mf{K.23') 


[(g'Aj')+i(gVj')]A:4 


(B.14) 
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4p = -5p5Ajp{(^Aj')/C3 + 

[{q'J^j') (B.15) 

4 a = 9pgAJp[iP'h')^ii - 

[{q'^j') + Kq'‘pj')]^3 \; (b.16) 

q± ) 

R% = gv9v\{PRj'){lC4j) + (PAj)(/Cij') - (jAj')(/CiP) 

+m/[(jAj')A:3 + {j(pj')JC4\ 

[{jKj') - i(j^j')][(/CiP) - mfJC,][{qRq') + -^qg^q')] (b.i7) 

^_L^_L 

[JCij) [{q'Rj') - i(g Vj')] - ^ (Aij') [(gAj) + i(gV5j)]| ; 
q^ qi J 

= 4y(9^-9'); (B.18) 

4y = 5y5A{(^4')(A2j) + (PKj)(lC2j') - {fkj'W^P) 

-m/[(jAj')/C4 + (jw')A3] 

+ 4^ [(jAj') - i(j(,7j')][(A:2R) + mfJC4\[{qRq') + i(g<pg')] (B-19) 

(/C2i) [(g'Aj') - i(g Vj')] - ^ (A2j') [(gAj) + i(gv3j)] j ; 

4a = 5y5A{(^'Aj')(A2j) + (P'Aj)(A2j') - (jAj')(/C2P') 
+™/[(jAj')A4 - (j'Pj')A3] 

+ 4^ [(jAj') +i(jVj')][(A2^’') -■^/A4][(9A9') 

^_L^_L 

+ ^ {K.2j) [{q'Rj') + i(g Vj')] + ^ (A 2 j') [(^Aj) - i(gv5j)]|; 
q± qi J 


(B.20) 



January 15, 2015 1:28 WSPC/INSTRUCTION FILE 

Kuzn’ Rum' Shlen' Gen' Compt' corr 


18 A. V. Kuznetsov, D.A. Rumyantsev, D. M. Shlenev 


= 5a5A{(^Aj')(/Cij) + (PAj')(/Cij') - {jAj')ilC,P) 


-mf[{jAj')IC 3 + {j(fj')IC 4 ] 

+ 4^ [0'Aj')-i0Vj')][(A:iP) + m//C3][(gA(3'')+i((3'(/?g')] (B.21) 

{JCij) [{g'Aj') - i(gVj')] - ^ [(^Aj) + j . 

Q± 9_l J 

Raa = Raa(<1 ^ -?'); ( B - 22 ) 

We note that the obtained results allow us to extract the limiting case n = 0. 
In particular, the amplitude Adyy containing the vector vertices only, coinsides 
after corresponding transformations with the amplitude of the Compton process 
in a strong magnetic field, calculated earlier in Ref. [1^] (see also Ref. [37] where the 
amplitude of the type M^y was considered for the case £'= £ = n = 0). In addition, 
it is easy to check that the resulting amplitudes containing the vector vertices, are 
manifestly gauge invariant. 


Appendix C. The set of expressions for the amplitudes for forward 
scattering 

For generalization of the results obtained in Ref. [5] to the case of magnetized plasma 
we consider the process of a coherent scattering of the generalized current j off the 
real fermions without change of their states (the “forward” scattering). We remind 
that in this case we mean under the generalized current j in the initial state only the 
field operator of a single particle, while the generalized current j' in the final state 
could be both the field operator of a single particle, and e.g. the neutrino current. 
In this case: £ = £', s = s', Kia = 2{pA)a, JC 2 a = 2{(pp)a, 

IC3 = 2Mi, IC4 = 0 . 

Since this is a coherent process, the total scattering amplitude is obtained by 
summing over all states of the medium fermions. We obtain the foilwing results for 
the summed generalized amplitudes: 

= (c-i) 

i,n—0 

( 2?(i) 2?(2) J 

^ J _ ^k'k _I_ ^kk' _ I 

\ip + q)fi -mj - 2j3n {p - q)f^ - m'j - 2/3n J 


where ff{E() = [l + exp(£’^ — p,f)/T] ^ is the fermion distribution function, T and 
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fjLf are the temperature and the chemical potential of plasma correspondingly, 

(C.2) 


75 ( 1 ) _ 


gsg'sjsj's {[(^Ap) + 2 / 3 ^ + 2TO7] 

= 41(9 ^ -q ); 


Ugp — — gsgpjsjp{q^p) In,e ^n-iy-l. j 


(C. 3 ) 


75(1) _ 

^vs — 


gsgvJsmf |[ 2 (pAj') + (qAj')] 


q -2 1 


(C. 4 ) 


[[(^Aj') + \{qipj')]In,lTn,l-l + [(<?Aj') - \{qq)j')]In-l/Tn-l,t-l] 

V 


~\J~^ [[(^Aj') + i{q(pj')]InJ,-l1n-lJ.-l + [{qAj') - i{q(pj')]In,t^n-l,e] 


77(2) _ 
^SV — 


gsgyjsmf |[2(pAj') - (gAj')] +XLiy_i] 


(C. 5 ) 


2B£ 

+ 1/ ~ [[(^Aj') - i{q<fij')]ln,eln,i-l + [(<?Aj') + i{q(pj ')]ln-l,eln-l,i-l] 

V ^_L 


+ W [[(qAj') - + [(gA/) -\-i{q(fj')]ln,£^n-l,£] i ; 


4 s = gsgAjsmf[2{pipj') + iqifij')] ; (C.6) 

41 = 41 ( 9 ^- 9 ); 

4p = ~gpgpjpjp |[(9Ap) + 2/3£] £ +21^_iy_i] 

-4/3Vn£I„,^J„_iy_i}; (C.7) 

4p = 41(9 ^ - 9 ); 

T>yp = = ~gpgv3pppf{q'P3 ) ; (C.8) 
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'^AP = -gpdAjpmf 


■- 7-2 _|_ - 7-2 

.n,i ' 


(C.9) 


+ [[(^Aj') + i{q(pj')]ln,lln,e-l + [{q^j') - i{qipj 

V q± 


+ Kqvj')]^n,e-l1n-l,i-l + [(gAj') - i{q(pj')]ln^tln-l,t] I ; 


ti(2) / . 

^PA = -9P9AJP'm-f 


{(^AjO 


t2 7"2 


(C.IO) 


I2I3£ 


[[(<?Aj') - i{q(pj')]ln,iln,e-l - [(^Aj') + i{qipj 


12/3n 


[[(^Aj') - i{qipj')]In,e-i1n-i,e-i + [(gAj') + i{q(fj')]ln,eln-i,e] >; 


T^^vl = 9v9v{ [ipAMPAj') + {PAj)ipAj') - (jAj')[2/3£ + (pAg)]J (C.ll) 
X + 4:(3Vnl{jAj')In,iPn-l,e-l 


1 2131 

n2 


{PAj)[{qAj') + \{qpj')] + {PAj')[{qKj) - \{qpj)\ 


1 2131 


(PAj)[(gAj') - iiqifj')] + iPAj')[{qAj) + iiqpj)] 


Pn-l,t-lPn-l,. 


I2l3n 

r,2 


{pAj)[{qAj') - iiqpj')] + {pAj')[{qAj) + i{q(pj)] In,eXn-i,l 


12^71 


{pAj)[{qAj') + iiqpj')] + {pAj')[{qAj) - iiqifj)] 


Pn-l,i-lPn,e-l 


+ [2I3£ + (pAq)] [[{jAj') + i{jipj')]I^j_-^ + [{jAj') - iijipj')]ll_^j] 
+ 4/3V^ [((?Aj)(gAj') - {q(fj){q(fj')]ln,e-iln-i,i}-, 


Pyv ~ Pyviq ~q^ 3 ^ j O ) 
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75 ( 1 ) _ 
^AV — 


9v9A{[{P^j){j'w) + (^Aj')OVp) 


(C.12) 


l2|3^ 


— [{Pvj)[{q^j') + '^{wj')] + {P‘^j')[{qM) - i{qvj)\]Pn,tXn,t-i 


^ [{P^j)[{q^i') - ')] + {Pvj')[{9^j) + i{qq^j)]]In-i,e.-iPn-i,i 

qi 


l2/3n 


[{pfMiq^j') - Kqpj')] + ipfj')KqM) + 


-\{{pfMiq^j') + Kqpj')] + ippj')[{qM) - Kqpj)]]Pn-i/-iXn,e-i 
y ^_L 

+{ppq) [[(jAj') + - [{jAj') - i(jv3j')]2:^-i,^]}; 


PvA = 9v9A{[iP'Aj)ij'pp) + iP'Aj')ij(pp) 


+ {jAj'){qpp) - mjijipj')] 


(C.13) 


^ [{P'pj')[iqAj) + iiqpj)] + iP'ipj)[iqAj') - i{qipj')]]I„^eXn,e-i 

qi 


+ 1 /^ [{P'pj')[i9Aj) - i{qipj)] + {P'ipj)[{qAj') + i{qipj')]]Tr,-ij,-iTn-i,(. 

v 


- 1 /^ [{PPj')[{qAj) - t{qpj)] + {ppj)[{qAj') + i{qipj')]]Tr,^eIn-i,(. 
V ^_L 




-{p<pq) [[(jAj') - i(j(/?j')] 2 :^ 7 _i - [(jAj') + }; 
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T^aa = 9a9a{ [{P^j){pM') + ip^j){P^j') - + ml + (pAg))] 


^ \Pn,i +^n-iy-i] + 4:P\/ni {jAj')Xn,ePn-i,e-i 


{PAj)[{qAj') + i{q(pj')] + iPAj')[{qAj) - i{q(pj)] 


{PAj)[{qAj') - i{qipj')] + {PAj')[{qAj) + i{q(pj)] 


(C.14) 



I2(M 


Pn,tPn,t—l 


Pn-\,l-\Pn-\,. 


/2/3n 

n’i 


(pAj)[(gAj') - i(gpj')] + (pAj')[(gAj) + \{qipj)\ Xn,iln-i,i 


i2f3n 


(pAj)[( 9 Aj') + i(9pj')] + (pAj')[(gAj) - \{qq}j)] 


+(M^2 +^71^ + {pAq)) [[{jAj') + i{j(pj')]Xlj,_^ + [{jAj ') - i{j(pj ')]Xl_^ f\ 
' [{qAj){qAj') - {qipj){qipj')]Xnj,-iXri-x\] 


91 

PfA=PAA{9^-q, 3^3')- 

We notice, that the expressions for amplitudes Aiys, Aivp, Mvv and Mav 
are manifestly gauge invariant. 
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